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Abstract 

We present both analytic and numerical results on the position of the 
partition function zeros on the complex magnetic field plane of the q = 2 
(Ising) and q = 3 states Potts model defined on Feynman diagrams (thin 
random graphs). Our analytic results are based on the ideas of destructive 
interference of coexisting phases and low temperature expansions. For the 
case of the Ising model an argument based on a symmetry of the saddle point 
equations leads us to a nonperturbative proof that the Yang-Lee zeros are 
located on the unit circle, although no circle theorem is known in this case 
of random graphs. For the q = 3 states Potts model our perturbative results 
indicate that the Yang-Lee zeros lie outside the unit circle. Both analytic 
results are confirmed by finite lattice numerical calculations. 
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1 Introduction 



The Ising model is known to be exactly solvable in one dimension and also in two 
dimensions in the absence of magnetic field. The magnetic field breaks the Z2 
symmetry which complicates a possible exact solution of the model. However, it is 
still possible to prove for arbitrary temperature an exact theorem ( see P]) about 
the location of the zeros of the partition function on the complex fugacity u = e^^ 
plane, where H stands for a complex magnetic field measured in units of = T. 
Henceforth, these zeros will be called Yang-Lee zeros. The theorem establishes that 
the Yang-Lee zeros are located on the unit circle = 1, k = 1,- ■ - n . Here 
n is the number of sites of the lattice whose details like coordination number and 
topology are not important for the proof. The symmetry under a change of sign of 
the magnetic field m — ^ 1/m is a key ingredient but it is not sufficient for the proof. 
The proof is a bit technical and strongly depends on the form of the polynomials 
Pn{u) in the partition function. A simpler proof would certainly be welcome. The 
authors of j2] have suggested an analytic way to find the location of the partition 
function zeros. It is assumed that for a system defined in a periodic volume V = L'^ 
with r different phases it is possible to write down the partition function of the 
system in terms of some functions (Z = 1, ■ ■ ■ , r ) as follows: 

Z = Y.qie-^''^' + O(e-^/^»e-^^0- (1) 
1=1 

where qi is the degeneracy of the corresponding phase / and fi is interpreted as its 
metastable free energy. The quantity Lq is of the order of the correlation length 
while / = min3fJe(/fc). Following [2], we can use a destructive interference condition 
in order to find the zeros of Z within 0{e~^^^''): 

3fJe = 3fJe f^^^^ < for all k^l^m (2) 

(3V{^mftl-^mr^^j) = ^ mod 27r (3) 

where ff^ = fi — {PV)~^^ogqi. It is assumed that the metastable free energies are 
functions of some complex parameter z which is in our case the magnetic field. 
The continuity of the real part of the free energy across the line of zeros already 
appeared in [3. It is clear that equations (j21) and © are useful whenever we have 
a closed expression for the free energies of the different phases of the system. In 
the case of the two dimensional Ising model in the presence of a magnetic field, 
although we do not have exact expressions for the free energies we can use low 
temperature expansions (LTE's) and the symmetry H — »• —H as in [2| to furnish 
a simple proof of the Lee- Yang theorem valid in the thermodynamic limit and for 
low temperatures. In this work we apply the same ideas for a particular case of 
fluctuating random lattice. For such lattices no circle theorem has been proven. 
The spins degrees of freedom are placed on vertices of Feynman diagrams which 
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are themselves dynamical degrees of freedom and need to be summed over in the 
partition function with appropriate combinatorial factors. In the case where the 
diagrams have double lines and are generated by a random N xN matrix model (see 
[4j^) their sum mimics, in the continuum limit, the Ising model on random surfaces 
(two dimensional gravity). Early EI and recent [H] numerical results on graphs 
of planar topology indicate that the Yang-Lee zeros lie on the unit circle although 
no circle theorem exists in this case. Notice that the form of the polynomials used 
in the proof of the Lee- Yang theorem is not preserved by their linear combinations. 
However, we believe that there might a generalization of the circle theorem for the 
new class of polynomials obtained by linear combinations with positive coefficients 
of partition functions of the Ising type. In fact, further support to this conjecture 
has been given in jS] where we studied the effect of the topology of the graphs on 
the position of the Yang-Lee zeros. We have noticed that changing from the planar 
to the torus topology the Yang-Lee zeros remained on the unit circle with a tiny 
change in their positions which is in agreement with the existence of an underlying 
robust theorem. Furthermore, in order to take the continuum limit we take large 
matrices N —>■ oo which allows one to find (see |3]) with help of the orthogonal 
polynomial technique a closed expression for the free energy even in the presence 
of the magnetic field. Thus, an analytic study of the free energy singularities is 
possible and leads to pure imaginary values for the magnetic field (see P), i.e., 
\uk\ = 1. Of course, this is only correct in the thermodynamic limit. On the other 
hand, the same issue of Yang-Lee zeros on a dynamical lattice can be studied in a 
simplified model of thin graphs where matrices become numbers ( = 1 ) and our 
final partition function becomes a linear combination of partition functions on single 
line Feynman diagrams. In this case no random surface interpretation exists and the 
model exhibits a ferromagnetic phase transition with mean field exponents, see ^U] 
and references therein. Once again numerical results for finite number of vertices n 
indicate the unit circle as the locus of Yang-Lee zeros P]. However, since we do not 
have a random surface interpretation, the results of [H] do not go through the case 
of thin graphs. In particular, the location of the Yang-Lee zeros is not known even 
in the thermodynamic limit. In this work we analyze this question starting from a 
closed expression for the free energy of the model obtained in the thermodynamic 
limit by means of a saddle point approximation both for the Ising and the Potts 
model. Then, we use the destructive interference equations of (2j to find out the 
location of the Yang-Lee zeros. Unfortunately the presence of the magnetic field 
complicates the form of the resulting expression for the free energy but an argument 
based on low temperature expansions and a symmetry {H —H) between different 
saddle point solutions will allow us to prove that the Yang-Lee zeros of the Ising 
model on thin graphs should be on the unit circle in the thermodynamic limit. The 
proof is exact for the whole range of temperatures for which the LTE's converge. For 
the case of the Potts model we have not been able to promote our low temperature 
results to an exact level but perturbatively we show that the Yang-Lee zeros should 
lie outside (but close to) the unit circle for low temperatures. The precise location 
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of the zeros seem to change with the temperature in a comphcated way. We also 
display numerical results in agreement with our analytic findings. 

2 The Partition Function on Thin Graphs 

We can define the q-state Potts model as a generalization of the Ising model (g = 2 ) 
where on each site of the lattice there is a spin degree of freedom a which can take 
q different values a = 1, 2, ■ ■ ■ , g. By summing over all spin configurations {cTj} on 
a given static lattice G„ of n vertices we obtain the partition function Zq{Gn) ■ In 
the presence of a static magnetic field H in the a = 1 direction we have, 

Z,{Gn) = E e^^<->'--"^ +'^^"-'^-' , (4) 

where Z]<jj> is a sum over all bonds ( propagators ) of Gn including loop bonds 
which connect a site ( vertex ) to itself. Here we are interested in the case where the 
lattice is a dynamical degree of freedom and the final partition function is obtained 
by summing over all Gn with n vertices: 

= EZ.iGl'^) (5) 

k=l 

In our case K(n) will be the total number of Feynman graphs G^^^ with n cubic 
vertices 0^ and no external legs. This clearly includes graphs of different topolo- 
gies. Each Zg{Gl^^) as well as Zj^'^^ will be proportional to a polynomium Pn{u) 
in the fugacity u = e^^ . We take cubic vertices since those are the simplest ones 
after the trivial quadratic case 0^ which corresponds to a D = 1 lattice. The parti- 
tion functions Zjj^'' can be generated by integrating over q zero dimensional fields 
Xi, - ■ ■ ,Xq representing the q different spin states a = 1, 2, ■ ■ ■ , g respectively, 



•7(2n) _ 



2m J g 



dg 

,2n+l 



I'D fie 



1 2 ^z>7 ^ J 3 ^ 1 Z_-^z — 2 * ^ 



V 



IV fie 



(6) 



where T>fi = Y\f^i dXi and T„ is a numerical factor independent of the magnetic field 
and temperature which is necessary to get rid of the bad asymptotic behavior of the 
undecorated 0'^ graphs. The contour integral is necessary to single out the term g'^^ 
containing 2n cubic vertices. Notice that vanishes for odd number of vertices, 
therefore we write the total number of vertices as 2n . The constant c will be related 
to the temperature. To each cubic vertex {Xf) representing a spin in the direction of 
the magnetic field there will be a m = e^^ factor in agreement with the paramagnetic 
interaction term in The quadratic terms in the argument of the exponentials 
in eq. (0) are responsible for the free propagators < XiXj > (or bonds < aiaj > ) 
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between the q different vertices associates with the states cxj. Taking ratios of the 
propagators and comparing with the ratios of the corresponding Boltzmann weights 
of (0)) we identify the parameter c with a function of the temperature as follows. 
The action in the numerator of eq. © can be written as 

So = \t X^K^^X, - f (e^^ Xf + ± Xf) , (7) 

where we have introduced the kinetic q x q matrix Kij = — c, ifi ^ j and Ku = 1. 
We can check that detK = (1 + c)''~^(l — {q — l)c). The propagators are given 
by the elements of the inverse matrix K~j^ = (det-ft')~^c ii i ^ j and i^j^^ = 
{det K)~^ [1 — (g — 2)c]. Therefore, using the Boltzmann weights from Q at van- 
ishing magnetic field we have 



e 



-/3 



{XjXj) _ c _ e^^i(^=°) 

JjQJQ ~ [1 - (g - 2)c] " e^-(^^=o) 

and consequently 

c=(e^ + g-2)-^ (9) 

Therefore corresponding to < T < oo we have respectively the compact range 
< c < l/(g- 1) . 

In the thermodynamic limit n ^ oo one can evaluate the free energy 

/r = (10) 

by means of a saddle point approximation as follows. The fist step is to rescale the 
zero dimensional fields Xi — > ^Xi such that Sg{Xi,u,c) {1/g'^) Sg=i{Xi,u,c) . 
Next, we decouple the contour integral over the coupling g from the integral over 



the zero dimensional fields by a change of variables g —>■ gy Sg=i . Finally, we have 



fq = lim 



n— »oo 



log / D/ie-"'°s5«=^ + O(-) 

2n J n 



:ii) 



Therefore the free energy per site is given in the thermodynamic limit by [TT] : 



U = \ log Sg=r (12) 

where Sg=i corresponds to Sg=i at a solution of the saddle point equations dxiSg=i = 
. Henceforth it is always assumed g = 1 and we will treat the q = 2 { Ising ) and 
q = 3 states Potts model separately. 
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3 Ising Model 



Aligning the magnetic field in the direction of ai, which is represented by the field 
Xi, and using the notation Xi = x and X2 = y we have respectively the saddle 
point equations dxS = = dyS : 



X 

y 



cy 
cx 



ux 

y' 



(13) 



Using those equations we can write down the action at the saddle point in a quadratic 
form: 

1 r 



S 



6 



2 I 2 
X +y 



2cxy 



From (|13|) we have 



uyil -yf + c{y- 1) + c 







(14) 

(15) 
(16) 



In the absence of magnetic field {u = 1) the solutions of (fTBjl simplify (see ^^1) and 
we have two low temperature solutions below the critical value Ccr = 1/3 and a high 
temperature solution valid for the entire range < c < 1 : 

1 



yLT± 

yHT 

Defining the magnetization 



m 



1 + cT J(l-3c)(l + c) 



ux 



[ux-^ + y-^ 



(17) 
(18) 

(19) 



we can check that ?/lt+ and yLT- correspond respectively to m > 1/2 and m < 1/2. 
Both solutions collapse at m{xHT,yHT) = 1/2 for c = 1/3. Substituting (fTTj) . (fTS|) 
and ()15|1 back in (|TH) we have: 



S 



LT 



6 



S 



'1-c) 



HT 



(21) 



Notice that both low temperature solutions have furnished the same action Slt ■ 

The presence of a nonvanishing magnetic field ( u 7^ 1 ) makes the solution of the 
cubic equation (fTB|) awkward and not very illuminating. We have solved (fTBjl instead 
as a Taylor expansion around c = (LTE). We found three expansions yLT± and 
yHT which reduce to (fTTj) and (fT5|) at u = 1. They furnish the following actions : 



S 



LT+ 



1 



c 

2^ 



3m3 



+ 



l)c^ 



2m4 



(22) 
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^LT- - g-^-^ ^ + 0{c) (23) 

~ 1+^2 c (1 + m2)c2 [u^-^u'' + iy 
^ HT — — — ~t~ ~t~ 

m 2^2 

(24, 

Clearly Slt+ and SLr- collapse into (pUj) at u = 1. For u > 1 {H > 0) fLT+ < fiT- 
and the system is magnetized in the direction of the magnetic field (x-direction ) 
while the situation is reversed for m < 1 (if < 0) as expected. In order to locate the 
Yang-Lee zeros of Z^^2 in the thermodynamic limit we first assume u = pe*^ and 
then solve equations (2) and (3): 



^e{fLT+-fLT-) = (25) 
3m hT^ = fLT+ + ^^^T^ (26) 

There are two basic ingredients which will allow us to locate exactly the Yang- 
Lee zeros for the whole convergence range of the LTE's without really worrying 
about the numerical details of the coefficients of the expansions ()22ll24p . First of 
all, notice that for i7 — — oo (or m ) the cubic equation (fTBj) has a unique 
solution y{u ^ 0) = 1 — and therefore there must be at least one LTE which is 
well behaved for u — > and terminate at the term proportional to c^. Of course, 
this can only be identified with Slt- which implies that their coefficients can only 
involve positive powers of u , i. e. , 

oo 

Slt-{u,c) = «n(M)c" . (27) 

ra=0 

Where each an{u) is a polynomial in u . The next important ingredient is a symme- 
try of the saddle point equations (fT^ which relate solutions with if > to solutions 
with H < . Namely, (|13|1 is invariant under 

u 1/u (28) 

Xi{l/u, c) — > uyj{u, c) (29) 

yi{l/u, c) uxj{u, c) (30) 

where {xi,yi) and {xj,yj) are, in principle, distinct solutions of (fT^. Back in (fT^ 
we are led to the relation: 

Si(^-,cj =u^Sj{u,c) (31) 
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The labels i,j indicate that the saddle point solutions on both sides of (jHH) do not 
need to be the same. From the first terms of the LTE's in (j22j) -(|24j ) and (|31|) we 
have the identification below 



Slt+{u,c) = \Slt-(-,c] (32) 



Sht{u,c) = ^Sht(-,c] (33) 



Putting back in the definition of the free energy and using (jTTj) we obtain 



2(/lt- - /lt+) = logM^ 

Coo OO 1 \ 

log J2 (^n{u)c^ - log E «n(-)c" (34) 

Using oo = 1/6 and u = pe^^, after expanding the logarithms we deduce: 

OO B{n) / , \ 

23fJe ihT- - fLT+) = 2 logp + ^ c" ^ Au,n cos{kd) - - (35) 

n=l fc=l \ P ' 

where B{n) and Ak^n are pure numerical factors. Thus, the condition ^e{fLT- — 
fLT+) = imply that the zeros must be on the unit circle Uk = e*^*= in the ther- 
modynamic limit. The condition on the imaginary parts of the free energies will 
locate the corresponding angles on the circle as a function of the temperature. 
We emphasize that our proof of the circle theorem for thin 0^ graphs hold inside 
the convergence region of low temperature expansions. 



4 q = 3 States Potts Model 

Using the notation Xi = x, X2 = y, X3 = z and fixing the magnetic field once again 
along the x-direction we derive the saddle point equations, dxiS = : 



X — c{y + z) = ux"^ (36) 
y-c{x + z) = y"^ (37) 
z-c{x + y) = z^ (38) 

From the difference of (j37p and (j38|) we deduce that there are two groups of solutions, 
either y = z or ?/ + z = l + c : 
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I) 



z = y (39) 
x = ^{l-c-y) (40) 
uy{y + c-lf + c{y + C-1) + 2c^ = (41) 



n) 

z = l + c-y (42) 
_ l±|l-4.c(l + c)]-/' 

- (1 + c)?/ + c(l + c + x) = (44) 

As in the Ising case it is instructive to first look at the solutions for vanishing 
magnetic field ( m = 1 ) . In this case we have from the first set of the solutions two 
low temperature and one high temperature solution : 



_ 1t[1-4c(1 + c)]V2 

yLT± — Zlt± — ^ 

.... ^ l + 2c±|l-4c(l + c)l./^ ^^^j 

xht = yHT = zht = 1 - 2c (46) 
Defining the magnetization 

Mg = 3) = ( 3^^"^' 3 ) (47) 
\ ux'^ + y-^ + z"^ / 

We identify the labels LT+ and LT— with solutions such that m > 1/3 and m < 1/3 
respectively. The solution LT— meets the high temperature solution with m = 1/3 
at c = 1/5. The second group of solutions PHj) does not furnish for m = 1 any new 
physical solution. Explicitly, for < c < (a/2 — l)/2, we have from the second 
group two real solutions: 



l-[l-4c(l + c)]i/2 
xi = yi = 7^ ^ (4J 
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1 + [1 - 4c(l + c)] 



1/2 



(49) 



and 



= ^^\^-Hi^c)Yr- ^^^^ 
y, = l + c+\[l-^c^l + c)Yl--c\ ^^^^ 

^ ^ l + c-||l-4c(l+c)p/^-c| (52, 

Notice that solution {x\^y\,z^ corresponds to a permutation ( ) of the 

solution LT+ given in ()45|) while the solution {x2,y2,Z2) coincides for < c < 1/5 
with the permutation {x,z) {z,x) of the solution LT—. However, its derivative 
is discontinuous at c = 1/5 and the solution deviates from LT— for 1/5 < c < 
(^2 - l)/2 . Concluding, we can safely neglect both solutions of the second group. 

Turning on the magnetic field {u ^ 1), the solutions of the cubic equation in 
()40|) become cumbersome. Once again we make a low temperature expansion around 
c = which we display below after the substitution in the action. From the first 
group of solutions {y = z) we have three possibilities: 

~ 1 + 2^2 2 + M c2 

^0 = ^r^-c^^ + ^(1 + 2^ + 3^2) 



3u 



^ - 13u^ - 12u^ + 2u + 2) + C(c^). (53) 



c 



2 



c3(2 + 3m) 



S, = i_,_,= + £!(H_M + 0(c') (55) 



Expanding the second group of solutions {y + z = 1 + c) we have two possibilities 

T 2 c3(m + 4) , 4, 



53 = g - - ^ + 0(c^) (56) 



l + M^ c c^(2 + u) 

u u 3 



-^4 = -^^^ + 0(c^) (57) 



Next we analyze the physical content of all those solutions. First, we notice that in 
the limit of vanishing magnetic field the only solution which correctly reproduces 
the LTE of Sht{u = 1;C) is 5*0 and therefore we will identify it with the LTE of 
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Sht{u,c) . However, the limit u ^ 1 can not be used to identify the remaining 
solutions since the solutions of the second group 5*3 and 5*4 become identical respec- 
tively to 5*1 and 5*2 for u — > 1 which on their turn correspond to the LTE of the 
physical solutions Slt+ and Slt- of the H = case ^. In order to understand the 
meaning of 5*1 — 5*4 we have to look at the limits H —>■ +oo ( u —>■ oo) and H —>■ — oo 
{ u ^ ). The only LTE solution well behaved for n ^ oo is 5*1 which will be 
therefore identified with a strongly magnetized low temperature solution, i.e., for 
this solution we expect limc^o'^ = 1- For m — >■ both 5*2 and 5*3 are well behaved 
and good candidates but 5*3 possess the lowest free energy ( recall / ~ log S ) as 
one can check numerically. Therefore 5*3 will represent the weakly magnetized phase 
which satisfies limc^o fn = . Thus, we have for the free energies of the physical low 
temperature solutions: 



2(/Lr- - /lt+) = log §3 - log Si 

= logM^+ (^^-2n-2^c^ (58) 

+ 3c^(^ -u^ + --2U-3] + 0{c') 
Vm u } 

We can find the position of the Yang-Lee zeros plugging u = pe*^ in ()58p and 
solving the destructive interference equations and (j^ . First, if we look at the 
leading terms ( zero temperature ) we will see from the equations and (f^H|) that 
the Yang-Lee zeros, in the thermodynamic limit, will be located on the unit circle 
p = 1 at T = 0. Our finite size numerical calculations are in agreement with those 
expectations (see Figure 1). Adding the next to leading terms of order we will 
still have the zeros on the unit circle but if we further truncate the LTEs at level 
the zeros will slightly move out the unit circle. Indeed, from (j25p and (j58|) we have 
at order : 

- (-) 

Imposing that — 1 < cos^^ < 1 we find numerically for < c < .5 that 1 < p < p* 
where p* increases with the temperature. It is worth commenting that the r.h.s. 
of is a monotonically decreasing function of p. Consequently, we have a closed 
curve outside the unit circle which is nonsymmetric across the imaginary axis with 
the farthest point from the origin being {p = p*,6 = tt). In Figure 2 both our 
numerical and analytic results are overlapped. The finite size results seem to tend 
to the (full line) analytic one in the thermodynamic limit. In Figure 2 we have used 
the truncation at order instead of although they differ by a tiny amount, 

^Incidentally, we notice that the discontinuous behavior of the solution H22|) only appears for 
1/5 < c < {\/2 — l)/2 which is outside the convergence region of the LTE of (|22|l . Remember that 
for < c < 1/5 both Slt- and (|22|l are identical, which explains why 5*4 and S2 become equal 
at u = 1 and low temperatures. 
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including the terms make the analytic results closer to the numerical ones. The 
terms are quite complicated to be displayed here explicitly. 




-1 -0.5 0.5 1 

Re (u) 



Fig. 1. Yang-Lee zeros of Z^J^^ at zero temperature, c = 0, for graphs with 2n = 20 (o), 
2n — 40 (x) and 2n — 200 (o) vertices. The closed curve is the unity circle. 



5 T ^ oo 

The results of previous two sections hold for the range of temperatures for which 
the LTE's converge. In general, we have not been able to go beyond those tem- 
peratures analytically. One exception is the case T ^ oo (or c l/(g — 1)) for 
which a subtle decoupling of the zero dimensional fields take place for arbitrary 
q similarly to what happens for the Ising model [q = 2), see jHI- The key point 
is to take a particular scaling for the coupling g appearing in the action ^ as 
c ^ l/(g — 1). Namely, the matrix Kij appearing in ((7j) has q—1 degenerate eigen- 
values Aj = 1 + c (j = 1, 2, ■ ■ ■ , g — 1), and the nondegenerate one Ag = 1 — (g — l)c. 
In order to decouple the fields Xj in the action it is natural to diagonalize the ma- 



12 



trix Kij through an orthogonal transformation^ Xi — > AijXj and rescale Xi — > 

The Jacobian cancel out in the calculation of Z^'^'^^ . After those changes the action 
becomes (repeated indices are summed over) 



q y2 

^ 2 

i=i ^ 



i. 
3 



u 



E 



AjkXk 



(60) 



0.5 



-0.5 



-1 - 




Fig. 2. Yang-Lee zeros of ^^^"g^ at c = 0.2, for graphs with 2n = 20 (o), 2n = 40 (x) and 2n = 200 (o) 
vertices, compared to the saddle point result (solid line) and to the unity circle (dashed line). 



Apparently, we have just transferred the mixing terms to the higher powers of the 
potential, however in the limit T — > cxo (c — >■ l/{q — 1)) we have ^ , while all 
other eigenvalues remain finite. Thus, after the redefinition g 
that limit : 



gXl^'^, we have in 



SAT 



2 



i=l 



(61) 



j=2 



^Explicitly the orthogonal matrix reads Aij = if i > j + 2 , Aig = and Aij = 

if i < i < g - 1 or j/^/j+f if i = j + 1 < q. 
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Thus, we have decoupled the zero dimensional fields on a specific point in the pa- 
rameter space while keeping information about nontrivial interacting terms. Clearly, 
this procedure is only useful if we do not care about the singular behavior of the 
original variable Xq. This is precisely the case of (0). Indeed, the Jacobian from Xi 
to Xi is singular in the limit T ^ oo but it cancels out because of the ratio in eq. 
((HI). Since Aig = l/^ we obtain 

SiT^oo) = ±^-^{y + q-l)Xl (62) 

Substituting this action in (jHI), the integrations over Xi, ■ ■ ■ , Xg_i cancel out leaving 
us with an expression which is finite in the limit c — > — 1) . From the term 
we have : 

")(T ^ oo) = T„(m + q - 1)2" (63) 

where T„ is a numerical factor independent of the magnetic field and temperature. 
Therefore we conclude that the Yang-Lee zeros of the g-state Potts model on thin 
graphs coalesce exactly at the point m = 1 — gasT— >cxo. The same result is valid 
on a static lattice ^Hl- Again our numerical calculations confirm this analytic proof. 

6 Conclusion 

We have proved that in the thermodynamic limit the zeros of the partition function 
of the Ising model on thin graphs lie exactly on the unit circle in the complex 
fugacity plane. Our proof is exact in the range of temperatures for which the low 
temperature expansions converge. For the case of the g = 3 Potts model we do 
not have the H —H symmetry anymore and our results were fairly perturbative 
in the temperature. In that case the zeros lie on closed curves which depend on 
the temperature. Those curves lie outside the unit circle and tend to the circle as 
T — i> 0. Our numerical results for small number of vertices seem to be in agreement 
with the analytic ones derived by means of the destructive interference formulas 
of Pj. The numerical results were very similar to the static lattice case treated 
in ^3] for both the Ising and g = 3 Potts model. We should mention that part 
of the motivation for this work came from a recent progress (see [TO] and jH]) on 
using the destructive interference formulas of j2] to find analytically the exact curves 
formed by the Fisher zeros (complex temperatures) of the q-state Potts model on 
thin graphs. However, the results of ITUj and [S] were derived for H = and the 
presence of the magnetic field complicates the form of the solutions such that we 
had to appeal to low temperature expansions. As a further work one might look 
at Yang-Lee zeros for complex temperatures as in |il6j as well as other types of 
vertices. Finally, we mention that, as in jU], we have also looked at connected 
partition functions obtained by first taking the logarithm of the ratio of integrals 
in © and then making the contour integral. The results for the position of the 



14 



Yang-Lee zeros were qualitatively similar to the ones reported here. We hope to 
return in the future to the open problem of proving the circle theorem for the Ising 
model on random lattices with finite number of vertices. 
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Figure captions 



• Fig. 1. Yang-Lcc zeros of ^^=3 at zero temperature, c = 0, for graphs with 
2n = 20 (o), 2n = 40 (x) and 2n = 200 (o) vertices. The closed curve is the 
unity circle. 

• Fig. 2. Yang-Lee zeros of Z^^^ at c = 0.2, for graphs with 2n = 20 (o), 
2n = 40 (x) and 2n — 200 (o) vertices, compared to the saddle point result 
(solid hue) and to the unity circle (dashed hne) . 
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